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SUMMARY 


This paper is concerned with the development of formulas of a recursive nature 
from a certain special type determinant. These formulas lead directly to computational 
routines that are effective in the evaluation of such determinants. Such computational 
routines are also applicable to a wide variety of recursion formulas, obtained by special- 
ization of the parameters involved. 

From the expansion of this type determinant, a recursion formula is obtained and 
later simplified to the form 
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This set of formulas is very efficient for computational purposes as discussed in the 
paper. 


I. INTRODUCTION 

This paper grew out of an investigation of an eigenvalue problem concerned with 
the Schrodinger wave equation using the Y ukawa potential. In that investigation, a 


recursion formula for the coefficients of the power series solution of the differential 
equation was used as a recursion formula and to deduce a determinant that determined 
the coefficients. It was discovered that routines for the effective evaluation of such types 
of determinants could be deduced as a generalization of that work and that these routines 
were more effective in precision and in time than the more general routines now avail- 
able. These routines may be used either for evaluations of these special type determi- 
nants or of quite general recursion relations. 


H. MATHEMATICAL TREATMENT 
The type determinant treated in this paper is the following 


Dn 


a ll 

a 12 

a 13 

• • • 

a in 

a 21 

a 22 

a 23 

• • • 

a 2n 

o 

a 32 

a 33 

• • • 

a 3n 



• 

• 

• 

o 

o 

o 


a nn 


in which a^, = O for i > j+1. By a straightforward expansion by minors, the following 
recursion formula is obtained. 
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where D 0 = 1. This formula will generate D n in n steps by taking in turn the values 1, 
2, . . . , n for i in equation 1. 

Formula 1 can be greatly simplified by setting (-1)* 1+ *Dj = a 21 . . . a. .R. for 
j=i, 2, . . . , n. This results in the formula 

i 


-a 


. , 4 .R. = y. a. . . .R. . 
1+1,1 1 j=i 1- J+M 1- 3 j 


(2) 


where R Q = 1 and for the last step in the process n = -1. Then it easily follows 

that 
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By changing the sign of a ^ . with the substitution b^ . = -a. + ^ . , we obtain the 
following set of formulas that are complete 
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It has been assumed in obtaining equations 2 and 4 that . 4= O for each i as is 
evidenced by the last equation in 4. 


If b. . . = O, then set 

i+l,i 


D - b 21 b 32 . . . b^ Yj 

j=i 


-j+l,i 


R. . 

i-3 


as is done when i = n. Proceed to calculate the remaining principal minor D' +1 as is 
done for D n . Then D n = D! +1 . This process may be repeated as many times as ze- 

roes occur immediately below the main diagonal. 


If zeroes occur (or are suspected to occur) , another procedure would be to base 
the calculation on equation 1, which is correct whether zeroes occur or not. In general, 
equation 1 leads to more computational time and effort than equations 4. 


It is evident from equation 1 and equations 4 that these formulas cover a wide 
variety of recursion relations and also are applicable in the solution of "almost diagonal" 
determinants. It is also true that these cases occur quite frequently in practice in 
physics, engineering, and other applied mathematics fields. These recursion relations 
might be termed "two dimensional. " For example, if we set a.. = a i+1 j +1 = a. for 1 > j 


and b , . = b., we obtain the 
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"one dimensional" recursion formula 
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which was associated with the Schrodinger - Yukawa differential equation mentioned 
above. 
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III. NUMERICAL CRITIQUE 


In approximating by evaluation the real roots of polynomials Pj, given by the 
determinant associated with relation 5 in the case just mentioned, the IBM 7090 required 
19 minutes for 645 evaluations of P 40 by the available routine. Since this technique was 
practically prohibitive, other procedures were sought. A procedure was devised based 
on relation 5 (actually programmed with the cooperation of Mr. Audie Anderson of the 
General Electric Company staff) . This procedure enabled the 7090 to make 645 evalu- 
ations (the same number) in 4. 2 minutes, not for just the case i = 40, but for all cases 
i = 1, 2, . .. , 75 (645 evaluations for each such case). This is convincing evidence of 
the importance of such procedures. 

It is of interest to compare the procedure outlined by equations 4 with the 
well-known process of replacing the a. + ^ . by zeroes or diagonalizing the determinant. 

To make clear the comparison, it is informative to write out a few steps for each method 
as is done below. 
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A close analysis of the two reveals several interesting facts. The two procedures 
are the reciprocals of one another in a certain evident sense. The procedure in the case 
of formula 4 calls for i multiplications and i-1 additions of known quantities in order to 
obtain Rp Here we count division as a multiplication by the reciprocal number as will 
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be done in both cases. The procedure in the case of the diagonal process calls for 
n i+9 multiolications and n-i+i additions of known quantities in order to obtain d x . At 
thTe„Znhe”du r e, in each case, there are n-i multiplications of *"own — es 
in order to obtain Dn- In the case of formula 4, obtaining R n calls for n-i multiplica- 
lions tstead of the expected n because = 1. On the other hand, the diagonal pro- 

cess only proceeds from i=2 to rtn. Taking all of these facts into account it is seen that 
precisely the same number of operations are involved in the two procedures, n^ y 
l(n-i) (n+4) multiplications and Jn(n-l) additions. On this basis, we have dea 

heat. " 

However, there is a difference between the two that might work against the dia- 
gonal process with respect to round-off error propagation. In the diagona procedure 
the divisors are rounded-off numbers; whereas, in the procedure of fornnda 4 ^he dm 
dends are rounded-off numbers while the divisors are given numbers Such a difference 

could be decisive in favor of the procedure of formula 4. Tests ar ®P^ e evaluate 1 
this difference on the IBM 7090 with determinants that are not too difficult t 

precisely. 


IV. CONCLUDING REMARKS 

From the simplicity of formulas 1 and 4, it is evident that programs can be 
written that are relatively simple because of the recursive properties ° f mSFC 

These programs are being written and will be made a part of the 1 rary 
Computation Division. 

Another point of interest and importance is illustrated by the difference in for- 
mulas i and 4. The two are mathematically equivalent as demonstrated above, but the 
difference in the number of multiplications is significant and, m general, formula 4 
would have a decided advantage in machine time. This illustrates the great importance 
of theoretical analyses of the machine routines that are to be use m prac ice. 
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